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The Asymptotic Representation of the Elliptic Cylinder 

Functions. 

By William Marshall. 



Introduction. 



We are indebted to E. Heine* for the first serious treatment of the functions 
of the elliptic cylinder. Starting from the defining equation 

« + (| COS 24> + 4 Z )^=0, 

where b and z are constants, Heine shows that E(q>) can be expressed in the 
form of the following infinite series : 

E($) = £a + <h cos 2$ + <h cos 4$ -{-.... , 

where the coefficients, functions of b and z, are subject to a certain recurrence 
formula, namely, 

«n + l = h ( n * — «)«n — a n-l • 

But this series, as was pointed out by Heine, is not convergent for all values of 
b and a, but only for certain particular values which he found as the roots of a 
certain function of b and z. This function, however, was not carefully defined 
by Heine, nor was his proof of the convergence of the series sufficiently rigorous 
to satisfy modern mathematical requirements. 

S. Dannacherf cleared up the inaccuracies and supplied the deficiencies in 
Heine's presentation. He subjected the function whose roots determine the 
convergence of the series to a careful investigation, and showed how these roots 
might be calculated. He gave also a satisfactory proof of the convergence of 
the series. 

* E. Heine, Kugelfunktionen, Bd. I, § 104. 
■f Inaugural Dissertation, Zurich, 1906. 
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W. H. Butts* made a particular study of the roots of this so-called limiting 
function, and calculated for different values of & a considerable number numerically. 

Iq the present paper, following an idea originally due to Stokes, f we obtain 
asymptotic or semi-convergent developments for the elliptic cylinder functions. 
These asymptotic expansions, which hold approximately for reasonably large 
values of the argument, have two decided advantages. In the first place these 
series are extremely well adapted to the calculation of the numerical values of 
the function when the argument is large, or even reasonably large; and in the 
second place they indicate without serious calculation the behavior of the function 
for large values of the argument, and particularly for an infinite argument, and 
they afford the most satisfactory method of determining where the function 
vanishes. 

I. The Transformation of the Equation. 

The equation of the functions of the elliptic cylinder may be written in 
the form 

< ^+(l?co S h z u + %)U=0, (1) 

where h and S3 are constants. If we change the independent variable by 

putting e u = z,| since 

, e u -{ e~ u 1 / , 1 \ 
cosh u = — = - ( z + - ) , 

the equation (1) goes over into 

d?U 1 . dU pY 2 l-S ,mrr- Q m 

We can remove the term containing the first derivative by putting Z7j = Z7z*. 
Then (2) assumes the form 



aru x 

dz % 



i+[S+(S+i+»)}+S>=o- « 



* Inaugural Dissertation, Zurich, 1908. 

t Trans. Cambridge Phil. Soc., Vol. IX, Part I, or Math, and Phys. Papers of G. G. Stokes, Vol. II, p. 3S9. 

% As a result of this transformation, or, in general, e aw = 2, the singular points of the resulting equation 
lie only at and 00 . Of this family of transformations, a = 1 is the only one which reduces the equation to 
a form such that, for large values of 2, it becomes approximately (9); and this form is desirable for our 
present purpose. 



the Elliptic Cylinder Functions. 



313 



Finally, if we put —=x, (3) becomes 



am 



dx % 
or, as we may write, 



+K?+i+*0^B-^=°. 



d%U * 4-/^1 4.P_|_?V7— 



(4) 



(5) 



where the following relations exist between the quantities involved in (5) and 
in (1): 



(6) 



p 




2 T 


i + »." 


2 


— 


& 4 
16' 




Vx 


=. 






X 


= 


2 * 





II. The Solutions in the Neighborhood of the Singularities. 
The point x = is an essential singular point of the equation(5), as is also 
the point x = oo . For if we put x = ^, (5) becomes 






(7) 



so that £ = is an essential singularity of (7). 

Following the known method, we might now obtain an asymptotic solution 
of (5) in the neighborhood of x = by putting 



Z7i = e x x a {a Q + a x x + a^t? + ), 



(8) 



and then determining a, /I and the coefficients a , a 1} a 2 , so that (5) is 

formally satisfied ; it will, however, in this case, be somewhat simpler first to 
obtain a solution in the neighborhood of x = oo , and from this to obtain the 
solution which holds in the neighborhood of x = 0. 
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For large values of x, (5) becomes approximately 

of which the complete solution is 

CTj = d cos x + <7 3 sin x. 



(9) 



(10) 



Our solution of (5), then, which is to hold in the neighborhood of x = oo must 
be such that for large values of x it assumes approximately the form (10). 
We assume, then, as an asymptotic solution in the neighborhood of as = oo , 

^ = sma(A + ^ + ^+....) + cosx(l? + ^ + |f + ....), (H) 

where the A 0f A 1} A % , ...., B , B u B z , are constants which must be 

determined so that the equation is formally satisfied. We have, then, from (11): 

U> = S mx(-^-^-^-....) + C osx(A + ^ + ^£+....) 
\ X* x 3 a; 4 / \ x xr / 

+ cos,(-|-^-^-....)- S in^ + f + § + ....). 

W = --.(^ + ^ + ?^. + ....) + »~..(-^-^-^— ") 
- 8in .( A+ ^ + A + ....) + eo 8 <i^ + ^ + ^ + ....) 

— in<-§-^-^-....)—, (*a§<§ + ....)- 

When we now substitute this and the value of Z7i from (11) in equation (5) 
and arrange according to descending powers of x, we have for the left-hand side: 



sin x 



1.2J-! , 2.3-A 2 , 3.4J. 3 , 4.5^4 , 
x s + a; 4 a; 6 "*" x 6 

. 2.1^ 2.25 a , 2.35 3 2.45 4 2^5A 



"^ a; 3 "^ x 3 "*" a 4 



+ -^- + 



a; 6 ' x 6 + 



o i M^-'-i . , ?A _ , 



X 5 + x 6 
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+ cos a; 



1.2-B, . 2.3.B, , 3.4B 3 . 4.5J3 4 . 

X 6 05* X 6 X b 

2.lA 1 2.2A 2 2.3A 3 2.4A i 2.5A 5 



+ 



A3 2 X S 

pB + pBj + 



X* 



X 



X 



X 



PB*- + P B * + 



X* 



+ 



X 

qB 1 



x 

pBt 

x s 



+ ... 



g^O . <l B l , 9 B 2 , 



X 



X" 



In order that this vanish, the coefficients of cos x and sin x must vanish 
separately. The necessary conditions for this are that the coefficients of the 
various powers of x, namely, x~ z , x~ 3 , ar 4 , . . . . , should vanish separately. 
This gives us the following relations for determining the A's and B's in terms 
of any two of them, say in terms of A and B Q : 



2 ° 

~- p A --AA, 



Ai = £ B = BJ x {p, q) (say), 

T» _ f A _ 

-t>i — — "2 ° — 



5,= 



2.2 

A(p + 2) _ 
2.2 



¥ ^p(p+2)=-A f 2 , 



B 

2 3 .2 



i>(p+2) = - J B / 2 , 



A -B f~ -ffGP + 2)(i> + 6) gl 
^ 3 ~^°L 2 3 7T! + 2.3_|~~ o/ " 

^ 3 ~ ° L 2 3 T3l 273 J ~~ A)/ » ' 

a - a r P (p + 2)(p + Q)(p+l2) g(j>+12) pq 1 _.. 

A *- A °l 2*741 2.4.2.3 ~ 2.1.2.4J- AoJi > 

^ ~ ^0 1_ 2*741 2.4.2.3 ~~ 2.1.2.4J ~ o/4 ' 

. B V p{p + 2){p + 6)(p + 12) ( j? + 20) (g + 12) (g + 20) 
5 °L 2 B .5! 2.4.2.3.2.5 

^~ — -°0/5 > 



.— pq(p + zo) _pq(p + z) 



(12) 



2.1.2.4.2.5 



2 5 .2 3 . 



B 5 = -A f 6 , 
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Here f u f 2 , f 3 , .... are simply abbreviations. We have, in general, 

A — B n~i{p+n{n— 1)] B n _ 3 q 
*»- 2(» — 1) ^ 2{n — lY 

K —- A n-llP + n ( n — 1 )l _ A-3? 

2(n — 1) 2(n — 1)'. 



(13) 



From these relations we may determine the constants A 1 , A % , A 3 , . . . . , 
B lt B 2 , B z , .... as far as we choose, though, as might have been expected, 
it is not possible to give general formulas for A n and B n except in terms of the 
coefficients which immediately precede them. Substituting now in (ll) the 
values of the coefficients as given by (12), we thus obtain for the asymptotic 
solution of (5) for large values of a; the following: 

L X X s X? X* J 

+ cos,[ J B -^_^ H -^ + ^-....]. (14) 

This we may write in the form 

\ X X Z X 3 J 

+ 5 (cos, + 4!i^_4^_^£ + ....). (16 ) 

If now we change the arbitrary constants by putting 

A sin x + B cos x == cos (a — x), j . ^ 

— ( A cos x — B sin x) = G sin (a — a;), J 

we may write (14) or (15) in the form 

n=0«.(.-.)[i-*+£-§+....] 

+ <7«n(a-<r)[^-^ + |-^+....] ; (17) 

or we may write 

U 1 = (7 [P cos (a — x) + 0sin(a — as)], (18) 

* This amounts to putting 

a = tan-' ^ • 
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where P and Q denote the following infinite series : 

(19) 



p -. J%_ I /t J§_ I 

x 2 "*" a; 4 a; 6 "*" ' 



H ~ x W ~ l ~ x 5 x 1 ~* ' ' " 

If we return to the original variables ?7and u, we may write (18) in the form 
(see the relations (6)) 

U=^ [Pcos(a-^) + Q sin ( a -^)], (20) 

where, in P and Q, x is replaced by -— and p and q by their equivalents. 

The series denoted by P and $ are at first, for reasonably large values of x, 
rapidly convergent; they are both, however, ultimately rapidly divergent, as 
may be seen from the relations (12). The point at which they begin to diverge 
depends on the relative magnitude of/ m and x n . In fact if we stop after taking 
In terms of the solution (14) (n terms containing sin x and n containing cos x) 
and set the result in (5), we obtain 

+ B co S (x + r ^7 t y]-, (21) 

and the right-hand side of this is small if "** is also small. In using such 

series in numerical computation, we need not take all the terms up to and 
including the smallest ; the number taken will depend upon the closeness of the 
approximation desired. 

To obtain a solution of equation (5) which holds in the neighborhood of the 
other singularity x = 0, we proceed most conveniently as follows : in (5) we put 

x = -j-, (5) becomes then 

t 

d % TJ l , 2 dZ7i , /l , p . \ _ n 

We can get rid of the term containing the first derivative by putting Z7i = U^t' 1 . 
Then (22) becomes 

^ + (J+|+ff)Di = 0. (23) 
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Finally, if we put t = i-*? - *, we obtain 



(24) 



Equation (24) is, however, of exactly the same form as (5), so that we may use 
the former solution (14) or (15). The solution of (24) which is valid in the 
neighborhood of £ = oo is then 



U z = sint(a + b f- a f- b f + .... 



V 



+ cos £ (b, 






(25) 



where the / a , / a , / 3 , .... have exactly the same meaning as in (12). If now 
in (25) we restore the values of the variables U z and £, we have 

a 6 /i a^ap/a __ x s h f s 



U 1 = x sin =*- ( o H ^ 

1 a; \ ° <r 



+ a;cos|^(V 



gQo/i _ ^0/2 , s 3 «o/ 3 



3" 



(26) 



We may write (26) in the form 

x 2 



TJ X — ancc-jsin — ( 1 



x 






5* <f £* 

+ Seisin |(-^-/ 1 -|-/, + |- /,-....) 

+~4(»-t* + f / «- ••••)[■ 



■)l 



(27) 



We may adopt then, as the two fundamental integrals of equation (5) which 
are valid asymptotically in the neighborhood of x = : 

U lp =x(Asm-^- — Bcos^), 

\ X X / 



w 



here 



U lq — x (B sin -2_ + A cos S—j f 



/v.2 ™4 ™6 

^=i~§- /.+ f-/.-fr/.+ 

t» *** ^ ^ / _i 83 /• 



(28) 
(29) 

(30) 
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These integrals have no meaning when q = 0. This case, however, we need not 
consider, since for q = (5) becomes 

tfJ' + (l+-£)Di = 0, (31) 

and this is a reduced Bessel equation with n z = — p.* 

It remains to bring the integrals (18) and (28) or (29) into connection ; 
that is, to determine the constants C and a in (18) so that (18) shall be an 
approximation to the fundamental integrals which are denoted by U lp and U lq 
in (28) and (29). In general this is brought about by expressing the function, 
here U lp or U lq} in the form of a definite integral and then, from this, deter- 
mining the leading term in the function as the argument increases. But the 
expressing of these functions U lp and U lq in the form of definite integrals of 
sufficient simplicity seems to involve difficulties not as yet overcome.f However, 
for any particular case, for any particular pair of values of p and q, the constants 
C and a might be determind in the following manner : Having chosen a value of 
x for which (18) and (28), for example, give a sufficiently close approximation, 
compute the numerical values of these expressions (18) and- (28) for this value 
of x. This would give one relation between G and a. By proceeding in a 
similar way for a second value of x, the values of G and a would be determined. J 

III. The Roots of U a — 0. 

Although we have not determined the constants Cand a in (18) so that (18) 
shall be an approximation to the fundamental integrals U lp and U lq , yet it is 
possible to determine from (18) the general behavior of these functions for large 
values of x. For the sake of simplicity we carry the discussion through with 
the variables TJ X and x and denote by U u either of the integrals U lp or U lq , 
where the a will serve to indicate that on the value of the constants a and G 

* That is, if in (31) we put U x = Jxk, we have, after reduction, 

1 



,„ + ^ +/ (l + ^_Jt) =0 ; 



and this is the ordinary form of the Bessel equation with re 2 replaced by —, p. 

tSee III, where the value of a is approximately determined, hut not in this general manner. 

JSee the numerical example in V. Here the arbitrary constant is determined in the above-mentioned way. 

42 



320 



Marshall: The Asymptotic Representation of 



will depend whether the one or the other of these integrals is intended. We 
have then, from (18), 

U u —G\_P cos (a — x) + Q sin (a — »)] ; (31) 

and we proceed to find where U la vanishes. We put 

P = M cos 4,, Q = Msm^, (32) 

where, as before, 

^~ X ~ x 2 "*" x* ~ x s + ••••' 

£) _ ./l__.ZL I jA_ /? .1 

v .c a; 3 ~ , ~ a; 5 a; 17 
(31) becomes then 

Ci,, = CJf [cos (a — a) cos ^ + sin (a — x) sin 4*], 
U u = CM cos (4> — a + x), 

where, from (32), since 

Jf = VP 3 + Q 2 , 



-♦^-i_S. 



tan 



4>-tan \— + ~ 3 + ^5 + 

/2 /"2 /" fi 

./2 9 /4 JlJS T 9 o 

Jf = 1 ^- 4- r^ - • • • 



a; 2 

Now, from (34), U la will be zero when 

cos (4> — a + ») = ; 
that is, when 



1 " 3" 

•>// — « + «: = -_-, -=- 



2n— 1 



7t; 



that is, when 



a; = a -J- 



2n— 1 



. 71 — 1^ 



(19) 



(33) 

(34) 



(35) 

.), (36) 

(37) 



(38) 



(where n is an integer). This equation we have to solve for x on the supposition 

2n— 1 



that x is large. If in (38) we put, for shortness, a + 

have to solve 

x = 4> — '»/'• 



71 = $, then we 



(39) 
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If now we expand •$> in terms of x by means of the expansion of tan -1 z, 
namely 

— z 3 4- — z 5 — — z 7 4- 



tan x z = z • 



we obtain, after reduction, 

/• / 1/2 /3 ~o~ /s /1/4 /1/2 + J2J3 / 1/2 H~ / 1/2 I c" 
+ = 4 + 3^ + 31 5- + ...- (40, 

Putting this value of $ in (39) and solving the resulting equation, namely 
* /1/2 f% 3" 

T X XT 

/s — fifi — /1/2 +72/3 /1/3 +/1/3 H £~ 
? - + ...-, (41) 

by the method of successive approximations, we obtain 

- = ♦-/,. i + (/,-/,/,+ 4 -/!). -£- 

(/s flfi /2/3 +/1/2 +/1./3 /H g~ 

+ 4/i/ 3 -4/?/ 2 + 4 + 2 /0^ + • * * * ( 42) 

We have, then, for the general formula for the rath root of the equation U la =0, 
the following : 

/3 

2 n 1 f fs A/2 H g - + / 1 

x = a + - ¥ -«- x + _ t -.... (43) 

a+—^-n (^a + ^-Ttj 

By means of the relations (12) we may express this in terms of p and q. We 
find, after some reduction, 

— 2n — 1 j) 5j> 2 + 6j>— 12g 

a= — a+ 2 ~ ^ — 2 a + (2n — l)rt + 3[2a 4- (2«— l)*] 3 ••■" W 

By means of the relations (6) this may be transformed into terms of the original 
quantities U, u, h and 33. 
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The expression (44) shows clearly how the roots are spaced along the axis. 
It shows also distinctly the influence of the undetermined constant on the 

position of the roots. The leading term in (44) is of course a -\ — n, so 

that a change in the value of a amounts to a displacement of the roots along the 
axis. (44) shows also how the position of the roots depends upon the values of 
p and q. In (44) q appears first in the fourth term, which for ordinary values 
of p and q, on account of the rapid convergence of the series, will be small, so 
that a change in the value of q will produce a very slight change in the position 
of the roots. 

In the calculations involved in the preceding pages — that is, in finding the 
values of 4 1 a s given in (36) and (40), and M as given in (37), and in solving 
the equation (41) — we may proceed as there indicated, performing the arith- 
metical operations upon the infinite series there involved, and carrying the 
results to the desired degree of accuracy. These results may be obtained more 
briefly (particularly if many terms are desired) and much more elegantly by the 
following method.* 

We had found, as a solution of (5), 

U la =C[Pcos(a — x) + #sin(a — »)]. (31) 

If now we put 

P = Mcos4>, Q = Msin^, (32) 

and further set C= 1, a = 0, we have, as a particular solution of (5), 

Ui = M cos ($ + x). (45) 

By substituting this integral in the differential equation (5), and determining 
the coefficients so that (5) is satisfied, we should be able to deduce recurrence 
formulas for the coefficients of M and ^. We have then, from (45), 

TJ[ — M' cos (as + $) — M sin (as + $) (1 + 4>'), (46) 

U{' = M" cos (a? + 4>) — 2W sin (a; + &) (1 + ty) 

— M sin (a; + $) W — Moos (x + 4<) (1 + ^') 2 - (47) 

*This method was suggested to me by Professor Burkhardt, who, however, has not published anything 
on the subject. It is generally applicable in similar problems, being particularly advantageous in the calcu- 
lation of the roots of Bessel functions. 
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When we substitute these results in (5), we have 

M" cos (as + 40 — 2if sin (» + 4)(1 + 4>') — if sin (as + 4) 4" 

- if cos (as + 4)(1 + V? + (l + f§ + Je) if cos (sb + 4) = 0. (48) 

In this the coefficients of sin (a; + 4*) and cos (x + 40 must vanish separately ; 
this gives the two equations for determining M and 4 : 

M<< — Jf(l +4/f + (l+K i + ^)M=0, (49) 

2 Jf (1 + 4') + M<V' = 0. (50) 

Equation (50) is readily integrated, giving 

1 + ^' = ^- ( & 1) 

When we substitute this value of 1 + 4' in (49), we have a differential equation 
for M, namely 

Jf "-jp + ( 1+ & + fO if=0 - (52) 

Now, since we are concerned only with a particular value of M, we may assign 
to the constant of integration c x any value we choose. As we shall see later, 
it will be most advantageous to take Ci = l. We know from the equations (19) 
and from the defining equation for M, namely M= VP 2 + Q*, that M must be 
an even function of x, and moreover that the first term is 1. We therefore 
assume 

* =i + !? + !? +? + ■••■ (53 > 

Then follow 

Jf/ = _?^._i^_«*!+ ...., (54) 

a; 3 ar x 1 x J 

M ,i — 3 -^ 2 + 5_ 1 4^ 7J>& 6 + / 55 ) 

a; 4 as" ar 
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Also, by the binomial theorem, 

^ x* r x" + a; 6 ^ ' 

where the a 2f a if a 6} are given by the relations : 



J_ 

M 3 



(56) 






3« 4 , 



n bk % a± — 7& 4 a a „-, 

— 6^a 6 — 8& 4 a 4 — 10& 6 a 2 

ag _ 



— o«Jr 



(57) 



If now we put the value of M" from (55), M from (51), and -=p from (56) 
in (52), we have, after putting c x = 1, 







a 4 




X B 




X S 


1 


«2 
X 2 


- + 


h 

X* 


- + 


A 8 

X 6 


- + 


X s 


- + 


p 

X 2 


- + 


Ph 

X* 


+ 


Ph 

x e 


+ 


x s 


+ 




+ 


q 

X* 


- + 


qh 

X 6 


- + 


qh 

X 8 


+ 


<h 




a 4 




«6 




a 8 





= 0. 



x" 



x* 



X 



(58) 



When the coefficients of the various powers of x are put separately equal to zero, 
the following recurrence formulas are obtained : 



4h = 

4& 4 = 
4& 6 = 
4& 8 = 


— P, 

-(p + 6)h~ 

— (p+ 20) & 4 

— (p + 42)k. 


-q + h, 

— qh % + b 6 , 

— qh + h, 


4&„ = 


— [p+(n — 


0(»-2)]&, 



j — qiCn-i "T" o n , 



(59) 



where b if b 6) b 8) .... are the same as a if a 6 , a S) 
only with the last term dropped ; that is, 
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. . as given in equation (57), 



b 2 =0, 



be,— - 


±hfh 
2 ' 

— 5^a 4 — 7& 4 a 3 


"6 ~~ 

bo — 


3 

— 6ko.a 6 — 8& 4 a 4 — 10& 6 a 2 


u 8 


4 



(60) 



We have then, in (59), a formula for any coefficient in the expansion of M 
expressed in terms of the preceding coefficients, so that any coefficient can be 
at once expressed in terms of p and q ; naturally, for the later coefficients the 
calculation becomes somewhat laborious, owing to the complicated character of 

a z, a i> a e, when expressed in terms of p and q. Below will be found a list 

of the first few coefficients calculated in this way. 



Ka 






P 

4' 

5f Sp 

32 ^ 8 



4' 
_ 15p 3 37j> 2 



28 



32 



195j? 4 , 611y 3 , 1821j? 



8 q \l6^J' 

315p 45p z q 



2048 



256 



128 



+ 



16 



128 



1 57 pq , 5q* 



IQhq 



32 



32 



(61) 



We may use these coefficients of the expansion of M to good advantage in 
the further computation, particularly in the task of solving by the method of 
approximation the equation 

sc = 4>— 4>. (39) 



To show how this may be done, we have from (51), after putting c = 1, 

J/ = J i 

Expanding M ~ 2 by the binomial theorem, we have 

J7 — C 2 -|- °i _i_ c « _i_ J2. _i_ 



(62) 



(63) 
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where c 2 , c it c 6 , 



are given by 
Co = — 2k 



2} 



2c 4 = — 3k«c z — 4& 4 , 

3c 6 = — 4&jC 4 — 5A4C3 — 6& 6 , 



(64) 



If we now integrate (63), putting the constant of integration equal to zero, 
we have 



v x Sx s 



5x 6 



' ' t 



(65) 



from which ^ can be expressed in terms of p and q if desired. In order to solve 
the equation (39), however, we may use the value of ^ as given in (65). We 
have then to solve 



* = ♦ + £ + A + A + A+... 



3a; 3 "■" 5a; 5 + 7a; 7 



(66) 



The successive approximations are readily found to be : 

»! = *, 







«* = * + £ + »_, 

A A 




-0- — °2 M z 

X. =A + °2 + d + 
4 ^ A 3 ' 


40.304 , c 6 
3 5 

A* 



(67) 



and by means of (64) and (61) these can be expressed in terms of p and q, 
giving the results previously obtained in (44). 

The fact that in (44) the value of the roots of U u = depends only on 
the undetermined constant a suggests the following method of determining 
approximately the value of this constant. We have 



U lp = x(Asm^l—B cos ^l\ 



(28) 
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where A and B have the values given in (30). If we put 




A = N cos 0, } 
B=Nsin$, J 


(68) 


(28) becomes 




U lp = xN (sin ^1 cos 6 — cos ^1 sin $) , 


(69) 


f7 ]p =a;JVsin(^S — fl) . 


(70) 


This will vanish if 




S in(^-0)=O; 


(71) 


that is, if 




"~fliL (n=0,l,2,8,....). 

+ nn 


(72) 



We find now, as in (40), 

*=f+f(/./-/.-4) 

+ -i(/i-/ 1 /.-/i/5+/i/i-/!/i+/W+ "§-) + ■■■■ 

We have now to solve approximately the equation 
*/q 



n7t + ~yT + "jT \hh—h 3-^ + • ■ • • 

on the supposition that a; is small. 

It will be more convenient to write (74) in the form 

r (\_ xfi , fix* x s ( f\ . fj % -h-f\ \ v \ 

Vq \nn~ q\nnf + q(nnf q* \ (nn) 1 + (nn) s )^ " " ) 

The successive approximations are : 



£C 



a5i = — - = I (say), 
1 nn u •"' 



-/ 1 



*8 = ^ / 2'( 
«3 



/l ^ 



nit 



{nnf J — 1 " 

(1 f 2/" s \ 

^n" ~~ "(nn) 5 " + ~^njr) = ?3 ' 



(73) 



(74) 



(75) 



(76) 



43 



328 



Marshall : The Asymptotic Representation of 



If we put this approximate value of x equal to the approximate value obtained 
in (44), we have an equation containing only a. This is 

p 5j? 2 + Qp — 12g 



, 2n— 1 
a H 5 — « = *» + 



2(a+^fli7t) 24 (a 



24 f a + 



2n — 1 \ 3 



«) 



+ 



(77) 



This gives, as successive approximations for a : 

In— 1 



«1 + 
<*3 + 



2 
2w— 1 



p 

7t==ln + ^L 



, 2rc— 1 __ 7 , P 

a 3 T o " — ™ ' 2? 



llp 2 + 6j9— 12g 



24Z* 



(78) 



This approximate value of a, which should hold for certain values of n neither 
too large nor too small, would indicate that the true value of a probably is a 
rather complicated function of the two parameters of the equation, p and q. 

In a similar way the value of a might be determined so that (44) would give 
approximately the position of the roots of the other fundamental integral U lq j. 

IV. The Boots of Ul = 0. 

In a similar way we may compute the roots of U[ a = ; that is, we may 
find the maxima and minima of U u . We have 



where 



^ = |{Pcos( a -^) + <2sin( a -^)], 



P= 1 

v ke u 



2 2 / 2 



T&e 
2% 






LH_4_ 



2% 

¥e eu 



+ ■ 



2 6 / 6 



k 3 e 3 



m 



ffll 



(20) 



(79) 



Now, since U' a exists, and indeed also in the form of an asymptotic expansion, 
we may differentiate (20) term by term, the results having of course the same 
meaning as the original expansion. We thus obtain 

rr> C\h u d • f te w \A ( *Ap/ £^V / he u \ 
UL ^^[-2 e u -PBm^a~—) + co S {a--)P >-\- cos (a--j 

+ sin(a-^) . ^]--2 f [pcos(a-^)+ Q.in(— *f)]. (80) 
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sin C a — <r)br + <2 -|J +C0S ( a -^A p — 2 -^r°- 

From (79) we have 
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p, _ 2 a .2/ 2 _ 2*.4/ 4 , 2*.6/ 6 



/feV 



&V< 



jfeV" 



• > 



(8 2) 
For convenience we may return to the former variable, putting, as before, 



qi 2/j , 2 s . 3/ 3 2 . 5/ 6 

v !«>« ~ iSrSm IAJm T • • • ■ 



/fee" 



We 



ifcW 



he u 



x = — , so that we have 
2 



x 3 ^ a; 4 " a 6 



P == 1 ZL 4- jL* A 4 



H /l JZ\fi ft I 

V — — — ^r -t -^-"-^r-t- 



£C 



£C" 



a* x* "*" a 6 x 8 ^ " " ' ' 



XX" x x 

Then the expressions in the square brackets in (81) become 
xP + g'-| = * -f -/—+ ....=i2(say), 



P. 

2 



P , -«<3-4r = -/i-i 



ar 



w 



cc 



We put now, for shortness, 



/1 + 2 = — 0i (i>, ?), ( or simply — ft), 

/3 + «•' 2 == ^ 3 » 

/i + 2^8 = &> 



= £. 



(19) 



(83) 



(84) 



(85) 
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Then the series in (84) become 

L 35 s £C 4 X* J' 



\_x x 6 x 6 x 1 J 



(86) 



The equation (81) becomes then, after dividing by x, 

Bx sin (a — x) ■+• S x cos (a — x), 



where 



R x = 1 



9% 4. _ 

a* ' '«•* 



i^ 
x 4 



^ 1 = lL_4 + 



a; 



a; 



06 

a; 6 



9% . 

a 7 1" • • • • 



(87) 



(88) 



As in III, we put 



JSj = Mx cos 4i , 

/Si = ifi sin ^ , 



(89) 



where M\ and ^i have exactly the same form as the M and ^ in (37) and (40), 
except that the (/-functions now replace the /-functions. Then (87) becomes 



or 



M\ [sin (a — as) cos $1 + cos (a — x) sin 4>i] = 0, 
sin (a — x + ^1) = 0, 



which vanishes when 

a — x + ^ = — nn (n an integer or 0) ; 

x = a + nn + i^. 



that is, when 



(90) 
(91) 

(92) 
(93) 



We may solve this equation in exactly the same way as in III we solved 
the equation (39); we need only replace each /-function by the corresponding 
(/-function, due regard being paid to the algebraic signs. We thus obtain 



x = a -4- nn + 



9i 

a + nn 



9» — 9i9a + & +9i 
(a + nn) 3 



+ ••• 



(94) 
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If we express this in terms of the /-functions, we have 

1 



x = a + nit — 



a + nn 
or, in terms of p and q, 



-J- ; ; -^5 r 



(a + nnf 



(95) 



x — a + nn ■ 



__2 24 _ 24 ^ H 

a + nn (a + nn) 3 



5 4 . 103 3 
-48^+960^ 



+ 



If we put for x its value x 



169 
240 



3 107 113 , /23 , 11\ 

^--2T^-480 + <20^ + Tj 



2 



(a + «7t) 6 
, we have, finally, 



(96) 



1 

¥1 



a + nn- 



pJ-J 

2 
a + nn 



24* +24^ + 9 
(a + nnf 



+ 



(97) 



This formula gives, then, the position of the (w + l)th maxima of the curve U a . 
Here also a change in the value of the undetermined constant a has practically 
the effect of displacing the maxima or minima along the axis. A comparison of 
the equations (96) and (44) shows further that the maxima and minima as given 
by (96) alternate with the crossings of the axis as given by (44). 

The magnitude of the maxima and minima ordinates, or indeed of any 
ordinate, is of course given by 



u - =1 [ p «- (« -t) + « sin (— t)] - 



which may be written in the form 
where M is given by (61) as 



= 1 



P 



q*e 



Su 



+ 



p z + 6p — 4q 



and 4 1 is given by 



qe 



iu 



+ 



(20) 



(98) 



(99) 



4> = tan -1 



/ £ 



p 



-£r-2p + Sq 



qi&to 



+ 



13 6 5 i 
6</-12* 



10 



+ ■•••/ 



(100) 
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The length of any ordinate is then given by (98), thus depending upon both the 
undetermined constants and a. This formula also shows how, for increasing u, 
the lengths of the maximum and minimum ordinates continually decrease. 

V. Numerical Verification. 
The equation considered by Heine, Dannacher and Butts is written in the 



form 



d 2 E . / 8 



dip 2 



f (JL cos 2$ + 4z)E = 0. 



If we put $ = iu, we have 

d 2 E / 8 



du 2 



(A cosh 2m + 4z)E = 0. 



(101) 



(102) 



If now we make the following substitutions : 

cosh 2u = 2 cosh 2 w — 1, 



16 



= ¥, 



(103) 



equation (102) becomes 



*B+4- = 4*, 



^- (^ cosh 3 « + $)#= 0, 



(104) 



which differs from equation (1) only in the sign of the second term. Now, by 
means of transformations similar to those employed in I, we reduce the equation 
(104) to the following form : 

d 2 E, 



du? 



-*--?■)*=•. 



where 



x 2 x i 

J ¥_ 

4 2 



», 






16 



(105) 



(106) 



For large values of x, (105) becomes approximately 

-drf—-® 1 - ' 



(107) 
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so that we may assume as an asymptotic solution of (105), for large values of a, 

E 1 = e°(A + ^ + ^- + ^- + ----} (108) 

Proceeding now as in I, we find, for the semi-convergent solution of (102) for 
large values of x, 



E—A, 



he u 


[~i l hl I hz I *" I 




11 Jce u ' /ke u \ 2 ' //^e u Y 

2 V 2 y \ 2 y 



(109) 



where Aq is an undetermined constant, and where the h lf h % , A 3 , are abbre- 
viations for the following : 

A — r 



_ r(r + 2) 

2 ~ 2 2 . 2 ! ' 



7i s 



r(r + 2)(r + 6) _s 



2 3 . 3 ! 



2 . 3 



* _ K-ijr + n(n — 1)] + ^-s 
rt »— 2(n — 1) 



(110)* 



and where 






A 3 



4 
16 * 



•8, 



(111)" 



Also between the constants of equation (102) and (109) there exist the following 
relations : 



r = — -£ 4a, 

16 



— ~ yi , 

4 
V6' 



& = 



(112) 



*It may be remarked that \, \, \, .... are the same as f lt f it f a , ..... as defined by (12), except that 
here all the signs are positive. Also p — — r and q = — s. 
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A solution of (102) as given by Heine is 



E{u) = -—- + a t cosh 2m + <h cosh 4u + , 



(113) 



where a = 1, a x = -5- 6a, and the other a's are given by the recurrence formula 



a n+1 = b(n 2 — z)a m — a Vr . 1 . 



(114) 



This series (113) converges (for any assumed value of b) only for values 
of z which are the roots of a certain transcendental function.* For large values 
of u it at first diverges rapidly, and ultimately converges slowly. The series 
given by (109), on the other hand, at first converges rapidly, for reasonably 
large values of u, and then diverges rapidly. It may be used for approximate 
numerical calculation, provided care is taken not to include any terms after 
divergence begins. We wish to show by means of a numerical example how 
the asymptotic expansion given in (109) leads to the same result as (113), but 
with much less labor in computation. 

It has been shown by Butts "f that if b be taken equal to .1, then one root 
of the limiting function, that is, one value of z for which the series (113) con- 
verges, is 5.58134 (correct to five decimal places). He gives the values of the 
coefficients computed for this value of b and z as 



«<> =n 


1.0000, 1 


Ox =- 


- 0.2790, 


«2 = — 


- 0.8721, 


«3 = 


0.4169, 


a 4 = 


1.0147, 


«6 = 


0.6402, 


«6 = 


0.2284, 


«ij = 


0.0548, 


«8 = 


0.0096, 


«9 = 


0.0015, 


«10 = 


0.0001. 



(l")t 



*See Dannaoher, I. c, p. 29. 

+ W. H. Butts., I. «., p20. 

%a„ and a 10 are not exactly as given by Butts, but have been changed slightly by means of Heine's 



formula a. 



'n+l 



< 



6n* 



which holds from a certain n on. 
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If we put now u = .7, equation (113) gives us 

E(.7) = .5 —.2790 cosh 1.4— .8721 cosh 2.8 + .4169 cosh 4.2 
+ 1.0147 cosh 5.6 + .6402 cosh 7.0 + .2284 cosh 8.4 
+ .0547 cosh 9.8 + .0096 cosh 11.2 + .0015 cosh 12.6 
+ .0001 cosh 14. (116) 

With a four-place logarithm table this gives 

E{.1) = .5000 — .6000 — 7.7980 + 13.9100 + 136.9000 + 351.0000 

+ 507.3000 + 494.2000 + 354.2000 + 222.5000 + 60.7500 = 2133. (117) 

In order to use equation (109) we must in (112) substitute the values of b and z, 
namely .1 and 5.58134. This gives 



r = — 22.47536, 
s = — 1600, 
h = W~io. 



(118) 



Substituting these in (109), putting «= .7, and using four-place logarithms as 
before, we have, using four terms of the series, 

E(.1) = 239500 A [l — .8824 + .3546 — (.1291 + .0765) 

+ (.0079 + .0138 + .0853)] = 239500^ (.3727). (119) 

A comparison of (117) and (119) serves now to determine A . We have, using 

four-place logarithms as before, 

^ = .02388. (120) 

This value of A may now be employed in finding the value of E(u) for other 
arguments, by means of the series (109). For example, if we put u = .S and 
use, as before, four terms of the series and four-figure logarithms, we have 

.#(.8) = 869600 A [l — .7976 + .2904— (.0956 + .0566) 

+ (.0052 + .0089 + 0572)] = 869600(.02388)(.4119) = 8552. (121) 

We may verify this by means of (113). Carrying the series to ten terms and 
using, as before, four-place logarithms, we have 

E(.S) = .5000 — .8672 — 8.5130 + 25.3400 + 305.20 + 954.2 

+ 1766 + 2004 + 1599 + 1345 + 444.0 = 8434. (122) 

Considering the roughness of the approximation used, and particularly the 
possible incorrectness of a 6 , a 7 , a 8 , a 9 , a 10 , this may be regarded as a very close 
44 
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agreement. The fact that (122) is too small may be accounted for by the fact 
that the next few terms would contribute materially to the result. In using 
series (113) for large values of the argument it would be necessary to know the 
value of a and the corresponding values of a u a 2 ,. . . . with much greater exact- 
ness. This fact determined in a way the choice of the particular values u = .7 
and m = .8 in this illustration. For values of u much smaller than .7, (109) 
does not converge at all ; and for values of u much greater than .8, the use of 
(113) is too uncertain. That the example is a particularly unfavorable one is 
due to the fact that b was chosen so small and z large in comparison, which has 
the effect of making both r and s large. Under more favorable circumstances 
(109) would converge much more rapidly, and for smaller values of the 
argument.* 

* Since the above was written, my attention has been called to an article by R. C. Maclaurin in Vol. XVII, 
Pt. I, of the Transactions of the Cambridge Philosophical Society. In this article, which is entitled " On the 
Solutions of the Equation (A 2 + k*)ip = in Elliptic Coordinates and their Physical Applications," Mr. 
Maclaurin reduces the equation (A 2 + k')f = to the form 

{z * - X) -§- + ■ -£- + Wxi -p' ] y = °- 

He then obtains power-series solutions in the neighborhood of — 1, + 1 )0D , this latter being asymptotic. He 
suggests the possibility of using this asymptotic expansion in finding expressions for the roots of the function 
itself and of its first derivative, and calls attention to the importance of these roots in practical applications. 



